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We derive a general theory of effective properties in metasolids based on phononic crystals with 
low frequency resonances. We demonstrate that in general these structures need to be described by 
means of a frequency-dependent and non-local anisotropic mass density, stiffness tensor and a third- 
rank coupling tensor, which shows that they behave like a non-local Willis medium. The effect of 
non-locality and coupling tensor manifest themselves for some particular resonances whereas they 
become negligible for other resonances. Considering the example of a two-dimensional phononic 
crystal, consisting of triangular arrangements of cylindrical shells in an elastic matrix, we show that 
its mass density tensor is strongly resonant and anisotropic presenting both positive and negative 
divergent values, while becoming scalar in the quasi-static limit. Moreover, it is found that the 
negative value of transverse component of the mass density is induced by a dipolar resonance, while 
that of the vertical component is induced by a monopolar one. Finally, the dispersion relation 
obtained by the effective parameters of the crystal is compared with the band structure, showing a 
good agreement for the low-wave number region, although the non-local effects are important given 
the existence of some resonant values of the wave number. 


I. INTRODUCTION 

Metamaterials are artificial structures with unusual 
constitutive parameters not found in natural materialsji 
like negative compressibility;^ refractive indejs^Ti^ or 
anisotropic mass density.— These properties offer new in¬ 
sights for the propagation of classical waves, and a wide 
variety of effects and applications have been found, like 
cloaking shells;^^ super-lenses;^ optical and acoustical 
black holesiSrJ^ or gradient index lenses— 

Metamaterials for acoustic or elastic waves, also named 
metafluids or metasolids, respectively, have been mainly 
implemented by means of sonic and phononic crys¬ 
tals, which consist of periodic arrangement of inclusions 
in a fluid (sonic crystal) or elastic (phononic crystal) 
matrix—^ If the inclusion is properly chosen so that it 
presents low frequency resonances, these structures be¬ 
have like effective materials with resonant-like consti¬ 
tutive parameters which can be either positive, zero or 
negative!^ 

Phononic and sonic crystals are anisotropic structures 
in general, therefore they present anisotropic constitutive 
parameters. Then, it was demonstrated that metafluids 
present anisotropic mass density not only near a local res¬ 
onance, but also in the quasi-static limit although for 
the case of metasolids it has been assumed in general that 
the mass density is a scalar— Recently, some works 
have shown that elastic composites have to be described 
by means of the so-called “Willis form” of the constitutive 
parameters /^1^^ which include a tensorial mass density 
and an additional coupling tensor, and the dynamic ho¬ 
mogenization of phononic crystals has also shown that 
this general description applies to these structuresi^ 

In this work the low-frequency limit of phononic crys¬ 
tals is analyzed, and it is shown analytically that low- 


frequency resonances actually induce an anisotropic mass 
density, which however becomes scalar in the static limit. 
It is also shown that the usual assumption that the neg¬ 
ative mass density is induced by dipolar resonances is 
not necessarily true, in a similar way as was previously 
demonstrated for metamaterials for flexural waves by the 
authors in a recent publication,-^ Finally, the dispersion 
relation of the full phononic crystal is compared with that 
of a homogeneous material with the obtained effective 
parameters, and a good agreement is found in general, 
although it is also demonstrated that non-local parame¬ 
ters have to be considered. 

The paper is organized as follows: After this introduc¬ 
tion, Section |n] explains the homogenization method em¬ 
ployed here. Following, Section Hill explains how to apply 
perturbation theory to derive some important properties 
of metasolids in the low frequency limit. Finally, Sec¬ 
tion |IV] describes a phononic crystal as a locally resonant 
metamaterial and Section |V] shows a numerical example 
of application of the theory. Section IVII summarizes the 
work. 


II. HOMOGENIZATION OF THE PERIODIC 
MEDIUM FROM THE BAND STRUCTURE 

The equation of motion of an inhomogeneous solid, as¬ 
suming harmonic time dependence with frequency w, is 
given by the classical elastodynamic equation^ 

- p{r)uj'^Ui = djCijki{r)dkUi ( 1 ) 

being Ui the components of the displacement field and 
Cijki the components of the stiffness tensor. Hereafter 
we employ the summation convention in which two re¬ 
peated indexes implies summation over their possible val- 
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ues, to simplify notation. If he medium is homogeneous 
the dispersion relation is obtained by assuming plane- 
wave propagation with wavevector k = kn, with k being 
the wavenumber and n a unit vector parallel to de prop¬ 
agation direction. Under this assumption, the equation 
of motion becomes the well-known secular equation for 
elastic waves^ 

= k’^nuCijnjjUj, (2) 

with Ui being the components of the displacement field 
and Cij the components of the stiffness tensor in Voigt 
notation (see Ref. and the Appendix 1X1). The solution 
for the dispersion relation uj = uj{k,n) is therefore given 
by the roots of the determinant of the matrix T defined 
as 

fy = puj'^S^j - k'^nuCijnjj. (3) 

In a phononic crystal both p{r) and Cijki{v) are peri¬ 
odic functions of the spatial coordinates, then Bloch the¬ 
orem is applied and the Plane Wave Expansion method^^ 
can be used to obtain the dispersion relation uj = ui{k, n) 
as the solution of the following eigenvalue equation 

i^^PG-G'{uG')i = (k + G)iiCfy^ {k + G')jj{uG')j (4) 

where pcjGf^j and stands for the Fourier compo¬ 

nents of the mass density, stiffness tensor and displace¬ 
ment field, respectively and summation over repeated in¬ 
dexes is assumed including the Fourier indexes defined 
by the reciprocal lattice vector G. The matrix elements 
(k + G)ii are defined in Appendix lAl In matrix form the 
above eigenvalue equation is expressed as 

uj'^Ngg'Ug'= Mgg'Ug'- (5) 

where 

{NGG')ij = PG-G'5ij, ( 6 ) 

{MGG'h = (fc + G)uCfj-^\k + G')j,. (7) 

This equation solves for the dispersion relation inside 
the phononic crystal, however in its current form it is 
difficult to figure out any property of the crystal as a 
composite. The description of the crystal as a material 
can be obtained by averaging the components of the dis¬ 
placement vector in the unit cell and finding in this way 
an equation similar to equation in which the coeffi¬ 
cients of the different terms multiplying the wavevector 
and the frequency define the effective parameters. The 
average of the displacement vector is given by the G = 0 
component of UG^ so that it can be obtained by express¬ 
ing equation ([SJ as 

lo'^NqqUq -\- co'^Nog'Ug' = MqoUo + Mqg'Ug' (8a) 
+ ^^Ngg'Ug' = MgqUo + Mgg'Ug'- (8b) 

It must be recalled that repeated indexes means summa¬ 
tion over all their possible values, and that hereafter it 


is considered that matrix elements labelled with G does 
not include the term G = 0, which is extracted from the 
above decomposition. We can now solve from the second 
equation for ug, 

ug' = —{Mg'g — ^'^Ng'g)~^{Mgo — w^Vgo)mo (9) 

and insert it into the first one, obtaining the following 
equation 

[lo'^Nqo — uj‘^Nog'Xg'g{Mgo — w^A^go) 

—Moo + Mog'Xg'g{Mgo — w^Vgo)] mq = 0 (10) 

where we have defined 

Xfry^i^-, fc) = {Mg'G - W^Vg'g)^™- (11) 

Equation m is formally the same as equation ®, 
however it is not an eigenvalue equation, but a secular 
equation for Uq similar to equation ([2]), where the so¬ 
lutions UJ = uj{k,n) are obtained from the zeros of the 
determinant of the matrix T defined as 

r = uj'^Noo — uj'^Nog'Xg'g{Mgo — w^-^go) 

— Mqo-I-Mog'XG'g(Mgo — w^Vgo)- (12) 

The matrix T is actually a 3 x 3 matrix, in which co¬ 
efficients are in general functions of both uj and k, what 
makes it less suitable for band structure calculation than 
equation ([S]) but more suitable for the description of the 
phononic crystal as a composite. Effectively, we can see 
that the elements of the Noo,Nog' and Ngo does not 
depend explicitly on the wavevector fe, 

(IVoo)ij = pSij, (13) 

(A^OG')ij = P-G'Sij, (14) 

(IVgo)^ = PGSij, (15) 

while the Mqo, Mqg' and Mgo contains this dependence, 

(Moo)y = hiCijkjj, (16) 

(MoG'h = kuCfj-^' (k + G')jj, (17) 

{MGQ)ij = {k + G)iiCfj ^ kjj. (18) 

The dependence with the wavevector and frequency can 
be reorganized then the T matrix can be cast as 

Tij = uj'^p*j - k'^miGijnjj - ujk{nijSij + SljTijj) (19) 

where the coefficients p*j , Cj j and Sjj are given by 

fe) = P^ij + ^‘^P-G'X?j^{^, k)pG (20a) 

C*jj{uj,k) = Cij- 

Gj^'{k + G')LiX?^^{uj, k){k + GUmG^j 

(20b) 

Sij(uj,k) = ujGj^'{k + G')uxfj^{^.k)pG (20c) 
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Equation (fT^ is similar to equation ([2]), but the con¬ 
stitutive parameters required to describe the phononic 
solid are more complex. This equation shows that the 
phononic crystal is a non-local Willis medium; ^^'^^ in 
which the mass density is a tensorial quantity and with 
the presence of the coupling field Sjj. The above ex¬ 
pressions are valid at any frequency and wavenumber, 
however in this work we are specially interested in the 
low frequency limit, that is, the limit in which the wave¬ 
length of the field in the background is larger than the 
typical periodicity of the crystal and it is described as 
a homogeneous material. It will be shown that even in 
the low-frequency limit these systems can have some spe¬ 
cial resonances in which the crystal behaves like a Willis 
medium with resonant and non-local parameters. 

If in equations (001) the limit a; —>• 0 and fc —>• 0 is taken 
it is found that the coupling field Sjj = 0, since it is di¬ 
rectly proportional to w. Also, the mass density becomes 
a scalar which is simply the volume average pij = pSij , as 
it is well known from the theory of composites. Finally, 
the effective stiffness tensor is given by 

Ch = ClJ - (21) 

and the medium behaves like an effective homogeneous 
medium with local and frequency-independent parame¬ 
ters. The mass density is a scalar, therefore all the in¬ 
formation about the microstructure of the composite is 
contained in the Gjj tensor, whose symmetry will depend 
on the background, inclusions and lattice symmetry. In 
the above expression it has been assumed the limit in 
which the frequency and the wavenumber tend to zero, 
however in practice this limit will be valid from zero to 
some cut-off frequency in which it will not be possible 
to neglect some terms containing the frequency or the 
wavenumber, the medium begins then to be dispersive 
and the constitutive parameters will depend on both fre¬ 
quency and wavenumber. It can happen however that 
these parameters be frequency-dependent even in the low 
frequency limit, under the condition what is called a “lo¬ 
cal resonance”. This happens when the parameter xgg' 
is singular, then it is found that all the constitutive pa¬ 
rameters can become resonant and locally singular, with 
the remarkable result that it presents, depending on the 
lattice symmetry, anisotropic mass density. 


III. PERTURBATION THEORY IN THE LOW 
EREQUENCY LIMIT 


we call resonances because then the x matrix is singular. 
Let us try to understand the nature of these resonances. 

Let us assume that we know the eigenvalues A„ and 
eigenvectors of the matrix M — lo'^N. We know then 
that the reciprocal of this matrix can be expanded by 
means of the eigen-decomposition theorem, thus we have 
that 


X={M-uj^N)-^ 


E 


vj, 8 ) Vn 
An 


(23) 


given that M — uS^N is actually a Hermitian matrix. The 
matrix M is defined in equation 0, and it can be ex¬ 
pressed as 

M = Mo + kMi + (24) 

where 

Mo = GuGfj-^'G'j^, (25) 

Ml = nuGfj-^'Cjj + GuGfj-'^'njj, (26) 

M2=nuGfj-^'njj. (27) 


For low frequencies and wavenumbers, the matrix M can 
be considered a perturbation of the Mq matrix, so that we 
can apply perturbation theory to relate the eigenvalues 
Afe with frequency. Let us define Un and G^jc? (with 
a being a quantity with units of length, for convenience 
in the units) the eigenvectors and eigenvalues of the Mq 
matrix, respectively, thus 

= G^Ja'^Un, (28) 


if we assume that fcMi -l-fc^Ma —uPN is a perturbation of 
the matrix Mq, the eigenvalues and eigenvectors A„ and 
Vn will be given, up to first order in perturbation theory, 

by 

A„ = C°/a2 + (29) 

Vn = Un + k'^ Ul + k^'^ anlUG 

i i e 

(30) 


where (assuming Un ■ = 1) 

Pn — 

G^n'>/a = u\MiUn, 
G^'> = ulM2Un, 


(31) 

(32) 

(33) 


The origin of the resonant parameters of these struc¬ 
tures is the X matrix defined in general as 

X = iM-uj^N)-\ ( 22 ) 

It must be pointed out that, in the static limit this x 
matrix simply is the reciprocal of the matrix M, however 
the term ui'^N can make that the determinant of the ma¬ 
trix M — uj'^N be zero for some specific values of ui, which 


and, for n ^ £, 


_ u^Nug 
C0/a2-C0/a2’ 

(1) _ ulMiUi 

Glla^-G^la^' 

( 2 ) ^ ulM2Ui 

GHo? - Gl! 


(34) 

(35) 

(36) 
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which ensures as well that Vn ■ Vn = 1. 

The dependence in k of the eigenvalues An implies that 
the effective parameters will be non-local in general, how¬ 
ever, as it will be shown later, metamaterials are in gen¬ 
eral designed by means of “soft” scatterers, that is, it is 
required that the velocity of the waves inside the scat¬ 
terers be much smaller than that of the background, in 
other words, the components of the stiffness matrix are 
in general much smaller than the density. Therefore, as a 
first approximation, we can neglect the coefficients mul¬ 
tiplying the wavenumber and approximate A„ as 

Xn^Cljc? -U?Pn. (37) 


and allows expressing the y matrix as (neglecting the 
perturbative terms in u) 




«)G'i(Mn)G3- 

C0/a2 - 


(38) 


This interesting result shows that at the resonant fre¬ 
quencies = Cn/Pn the effective parameters become 
singular, and in the neighborhood of this frequency the 
can have positive, negative or zero values. These reso¬ 
nances are determined by the ratio of the eigenvalues Cn 
of the matrix Mq and by the perturbation term pn- Also, 
the coupling of these resonances with the different con¬ 
stitutive parameters is defined by the symmetry of the 
eigenvectors of the matrix Mq, as will be explained in 
the following section. 

It must be mentioned that the eigenvectors corre¬ 
spond to the eigenvectors of the matrix Mg, which ac¬ 
tually is the matrix Mqqi given by equation ([7]) but for 
fc = 0 and removing the terms corresponding to G = 0. 
These eigenvectors correspond to a physical system which 
is easy to understand: Imagine a phononic crystal in 
which the stiffness tensor is periodic while the mass den¬ 
sity is equal to that of the background. It is easy to 
see that the eigenvalue equation of this system at the T 
point, that is, for fc = 0 , is according to equations ([H), 


uj'^uq = 0 (39a) 

= AfQQ/itQ/. (39b) 


which for a; 7 ^ 0 has the only solution ug = ( 0 ,Mn), 
being the eigenvectors of the matrix Mq. This rela¬ 
tionship between the eigenvalues and eigenvectors of a 
physical system with those required to compute the ef¬ 
fective parameters suggest that other numerical methods 
more powerful than the Plane Wave Expansion method 
could be used to characterize these systems, after prop¬ 
erly Fourier transform solutions and elastic constants dis¬ 
tribution. Therefore, on the basis of the present theory, 
this work opens a door to a more efficient calculation of 
the effective parameters, whose discussion is beyond the 
objective of the present work. 


IV. EFFECTIVE PARAMETERS IN THE 
LOCAL APPROXIMATION 


The study of a local resonance is made in a regime in 
which the wavelength of the propagating field is larger 
than the typical periodicity of the composite. This hy¬ 
pothesis implies that in equations (1701) we can make the 
approximation k + G ^ G, so that in these equations 
the dependence on k disappears and the parameters, al¬ 
though frequency-dependent, are “local” and given by 


ptj (w) 

= pSij + co"^ P-G'xfrr^ h>)PG 

(40a) 

ch(.^) 

= GlJ — Cj^ G'l^ixfn^ h)GmMC^J 

(40b) 

Sijh) 

= coCjh Ghxfj ^h)PG 

(40c) 


where *^(^) is computed from equation (ED as 
X^^ioj,k = 0). A better insight in the properties of 
these parameters can be done by including the expan¬ 
sion of X given by equation (1381) . then we have 


+ '^uPp-G' 

n 

Chico) = Cij- 


iuhcijUnjcj 

Cha^-oo^Pn 


Pg, 


E ^-G'r</ i'^V)G'liUn)Gm ^ ^G 

CHa^-uohn 


C ! \ g-G' G' ('^n)G'l{Un)Gj 

' 5 l ( w )=> Glit^TToi - CT—PG- 


Cha^-oo^Pn 


(41a) 


(41b) 

(41c) 


We can now define the quantities 


{^n)i — hniGiPG ^ (^^) 

(Kin)/ — hn)GmGniMC^ j ^ (^3) 


and, given that any Fourier coefficient satisfies = 

Fq we get for the effective parameters the following ex¬ 
pressions 


) Ec 70 /a 2_^2 


Chiuo) = Cjj-Y, 


{BhijBhj 

C^ja? -oj'^pn’ 


Sijh) = 


(P^n)I i-^n) j 
C^jaP- - uP-pn 


(44a) 

(44b) 

(44c) 


The above equations relate the effective parameters with 
the properties of the eigenvectors and eigenvalues of the 
matrix Mq, as well as with their perturbations. Let us 
note that, for a symmetric unit cell, we will have that 
F(—G) = F(G), and it is easy to understand that for 
this type of systems we have that 


u-G = ±ug, (45) 

which implies two type of solutions, (A„)i 7 ^ 0 and 
iBn)i = 0 when it_G = ug and (A„)i = 0 and 
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FIG. 1. Phononic crystal studied in the present work. The 
system consist of a triangular arrangement of coated cylin¬ 
ders in an epoxy background. The cylinders consist of a lead 
core of radius and a rubber shell of radius Ts (see text for 
numerical values). 


{Bi)n 7 ^ 0 when = —uc- The former induces a res¬ 
onant mass density, while the latter induces a resonant 
stiffness tensor. The two cases implies that the Willis 
tensor Sjj is equal to zero, and this also implies that a 
symmetric system cannot excite simultaneously a reso¬ 
nance in the stiffness tensor and the mass density, that 
is, we cannot have double negative materials in this way, 
unless the different resonances be too close each other. 

A deeper insight into the properties of symmetry and 
non-symmetry of these resonances is beyond the objec¬ 
tive of the present work, in which we want to focus at¬ 
tention on the properties of the resonant mass density, 
however a future work concerning non-symmetric lattice 
will be prepared and published elsewhere. 


V. NUMERICAL EXAMPLE: RESONANT AND 
NON-LOCAL ANISOTROPIC MASS DENSITY 

Figure [T] shows the system to be studied in the present 
work. It consists in a periodic arrangement of coated 
cylinders in an epoxy background {pb = 1.18 Kg/dm^, 
Eb = 4.35 GPa and Vb = 0.37). The cylinders are made 
of a lead core {pa = 11.34 Kg/dm^, Ea = 16 GPa and 
Va = 0.44) of radius Va = 0.16a and a rubber shell 
{p, = 1.3 Kg/dm3, Es = 2.7E - 4 GPa and i^b = 0.499) 
of radius Vg = 0.4a, since this combination of soft-hard 
coatings is known to present low frequency resonances. 
The cylinders are arranged in a triangular lattice, in this 
way we expect the effective material to be transversely 
isotropic. 

Figure [5] shows the effective mass density tensor rel¬ 
ative to that of the background pb as a function of fre¬ 
quency as computed by using equation (I44al) . The sym¬ 
metry of the lattice makes that any second rank tensor 
will have only two components, one for the xy plane and 
another one for the 0 plane. It is seen how in the low fre¬ 
quency limit the two components are identical and equal 
to the normalized average mass density pjpb, as expected, 
however it can also be seen how they split as a func¬ 



FIG. 2. Effective mass density tensor for the proposed 
phononic crystal. Both the transversal component (blue con¬ 
tinuous line) and the z component (green dashed line) present 
a local resonance. In the vicinity of this resonance both com¬ 
ponents of the mass density become negative. 


tion of frequency and present two different resonances, 
so that the system behaves like an elastic medium with 
anisotropic mass density. Moreover, it can also be seen 
how these components are negative in different frequency 
regions. It is found that the effective stiffness tensor is 
nearly constant in frequency in this region, and that the 
coupling field Sjj is zero, as expected from the discussion 
in the previous section. 

Figure [3] shows the field distribution of the lower fre¬ 
quency resonances of the mass density tensor depicted 
in figure [2l the upper panel for the z component and 
the lower panel for the xy one, left panels show the real 
part of the mode while right panels show the absolute 
value. It is interesting to note that the xy mode has 
a dipolar symmetry, as it is commonly assumed in the 
literature^ while the z mode has monopolar symmetry. 
The fact that a monopolar symmetry could induce a neg¬ 
ative mass density behaviour was already found by the 
authors in a recent paper— in the study of flexural waves 
in thin plates. This result is consistent with the theory of 
elastic waves in plates, given that a plate with a periodic 
arrangement of inclusions is indeed a finite slide of the 
two-dimensional phononic crystal studied here, and this 
result suggest that the propagation of flexural waves is 
mainly dominated by the z component of the mass den¬ 
sity. This important result should be taken into account 
in the homogenization theory of plate metamaterials, al¬ 
though a deep insight into it is beyond the objective of 
the present work. 

Equations (gOl) show then that the phononic crystal 
can be described by means of locally resonant consti¬ 
tutive parameters, whose frequency dependence can be 
easily computed. The description of a phononic crystal 
as a frequency-dependent homogeneous material will not 
be valid for every wavenumber and frequency, and to de¬ 
termine these limits the dispersion relation obtained by 
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FIG. 3. Resonant modes inducing an effective negative mass 
density. Upper panels for uikal2-KCt = 0.0072, corresponding 
to a resonance in pz, lower panels for u}kal2TTCt = 0.0194, 
corresponding to a resonance in p± (see text for further dis¬ 
cussion). 


means of the constitutive parameters is compared with 
the band structure obtained from the eigenvalue equa¬ 
tion ®. Given that Sjj is zero for this example and 
Cjj is constant in frequency, along the FJf direction the 


dispersion relation for the effective material is 

UJ^P±{uj)Ux = (46) 

uj'^p±{uj)uy = klCeeUy, (47) 

uj'^Pz{oj)uz = (48) 

while along the FA direction, that is, along the z axis, the 
dispersion relation is (notice that in this case C 55 = CI 4 ) 

uj'^p^{u:)ux = klCl^^Ux, (49) 

uj'^p^{uj)uy = klCl^Uy, (50) 

uj^Pz{uj)uz = klC^^Uz- (51) 


FigurelU left panel, shows the dispersion relation along 
the FX direction (x axis) computed by means of the 
eigenvalue equation ([5]) (black lines) compared with the 
dispersion relation obtained by means of the constitutive 
parameters. Red and blue dots show the results for the 
xy modes, and it is seen that there is a good agreement 
between the eigenvalue equation and the effective mate¬ 
rial dispersion relation. The dispersion relation for the z 
mode (green crosses) is however different from the eigen¬ 
value equation and the effective material, and there is 
an agreement only for very low wavenumbers. As will be 
seen later, the reason for this disagreement is that the lo¬ 
cal description of the metamaterial is not accurate here, 
and it is required the inclusion of the non-local compo¬ 
nents, that is, the dependence on the wavenumber in the 
constitutive parameters. 

FigurelU right panel, shows similar results for propa¬ 
gation along the FA direction (z axis). It is shown here 


r X FA 




hxCL hzCi 

FIG. 4. Left panel: Dispersion relation of the phononic crystal 
along the FX direction (black lines) compared with those ob¬ 
tained from the effective local constitutive parameters (green 
crosses and blue and red dots). Right panel: Dispersion re¬ 
lation of the phononic crystal along the FA direction (black 
lines) compared with those obtained from the effective local 
constitutive parameters (green crosses and blue dots). 


that the xy modes, which are degenerate given that the 
crystal is transversely isotropic, are perfectly described 
by means of the effective material, however the z modes, 
corresponding to green crosses, agree only for very low 
wavenumbers. There are also a set of flat bands that 
can be fairly difficult to predict by means of the effective 
material parameters. The reason for that is that these 
modes occur only at a given frequency and correspond to 
very sharp modes, and although they are properly pre¬ 
dicted by the theory as a resonant frequency their 
effect is difficult to see in the constitutive parameters. 


The spatial dispersion of the z mode can be under¬ 
stood by means of the calculation of the non-local con¬ 
stitutive parameters using equations (OIH) . Figure[5]shows 
these parameters at a frequency wa/27rct = 0.0028, corre¬ 
sponding to a frequency in which the z component of the 
local mass density is negative. The upper panel shows 
the non-local pz as a function of the wavenumber along 
the FA and FA directions. It is clear that the origin of 
the non-locality is a wavenumber resonance, for which 
the major contribution will have its origin in the x 
trix. It is also seen that the C 33 component, responsible 
of the propagation of the mode along the FA direction, 
also becomes non-local, while the C 44 remains constant. 
Additionally, the S 53 and S 33 elements, which are zero for 
k = 0, appear as resonant components. The contribution 
of these spatial resonances is essentially to displace the 
opening of the band gaps, as can be seen from figure IH 
for which their influence is important before considering 
only the local theory. 
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FIG. 5. Non-local constitutive parameters related with the 
propagation of the z mode in the phononic crystal. Upper 
panel: z component of the mass density. Mid panel: C 44 and 
C 33 components of the stiffness tensor. Lowe panel: 553 and 
533 components of the coupling field (see text for details). 

VI. SUMMARY 

In summary, it has been analytical and numerically 
demonstrated that phononic crystals behave as elastic 
metasolids with anisotropic, resonant and non-local ef¬ 
fective parameters, with the remarkable result that the 
mass density is also anisotropic in general, although in 
the static limit this quantity recovers its scalar nature. 
Also, it has been demonstrated that the symmetry of 
the resonance inducing this behaviour is not necessar¬ 


ily dipolar, as it is commonly assumed, while it can also 
be monopolar. The non-local and anisotropic nature of 
the mass density has important implications specially for 
the study of plate metamaterials, since these structures 
are essentially finite slides of phononic crystals. It must 
be pointed out that the generality of the equations de¬ 
rived can be used to the homogenization of phononic 
crystals with more complex unit cells, with the objec¬ 
tive of achieving double negative metasolids. Finally, the 
theory can be extended to phononic crystals with piezo¬ 
electric inclusions, where resonant piezoelectric constants 
are expected. 
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Appendix A: Matrix Notation 

Through the paper Voigt notation for the indexes is 
used, in such a way that lower-case indexes run from 1 
to 3 and upper case indexes run from 1 to 6. Also, the 
wavevector is defined in terms of the Vn matrix defined 
as 

/I 4 0 0 0 14 Vy\ 

y = 0 Vy 0 y 0 y (ai) 

\ 0 0 V, Vy V, 0 J 

Therefore the matrix elements (k + G)ii are 


J 


f kx Gx 0 0 Q kz Gz ky -I- Gy\ 

fc + G = j 0 T 0 kz -\- Gz 0 kx Gx | (-'^2) 

y 0 0 kz -\- Gz ky Gy kx + Gx 0 J 


being therefore {k + G)jj the transpose of the above ma- normal vector n, 
trix. Similarly, the same matrix n,;r is defined for the 

( Ux 0 0 0 Hz ny\ 

0 Uy 0 Uz 0 Ux \ (A3) 

0 0 Uz riy Ux 0 j 
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